Epitranscriptome insights into Riccia fluitans L. (Marchantiophyta) aquatic transition using nanopore direct RNA sequencing.

BACKGROUND: Riccia fluitans, an amphibious liverwort, exhibits a fascinating adaptation mechanism to transition between terrestrial and aquatic environments. Utilizing nanopore direct RNA sequencing, we try to capture the complex epitranscriptomic changes undergone in response to land-water transition. RESULTS: A significant finding is the identification of 45 differentially expressed genes (DEGs), with a split of 33 downregulated in terrestrial forms and 12 upregulated in aquatic forms, indicating a robust transcriptional response to environmental changes. Analysis of N6-methyladenosine (m6A) modifications revealed 173 m6A sites in aquatic and only 27 sites in the terrestrial forms, indicating a significant increase in methylation in the former, which could facilitate rapid adaptation to changing environments. The aquatic form showed a global elongation bias in poly(A) tails, which is associated with increased mRNA stability and efficient translation, enhancing the plant's resilience to water stress. Significant differences in polyadenylation signals were observed between the two forms, with nine transcripts showing notable changes in tail length, suggesting an adaptive mechanism to modulate mRNA stability and translational efficiency in response to environmental conditions. This differential methylation and polyadenylation underline a sophisticated layer of post-transcriptional regulation, enabling Riccia fluitans to fine-tune gene expression in response to its living conditions. CONCLUSIONS: These insights into transcriptome dynamics offer a deeper understanding of plant adaptation strategies at the molecular level, contributing to the broader knowledge of plant biology and evolution. These findings underscore the sophisticated post-transcriptional regulatory strategies Riccia fluitans employs to navigate the challenges of aquatic versus terrestrial living, highlighting the plant's dynamic adaptation to environmental stresses and its utility as a model for studying adaptation mechanisms in amphibious plants.

Assessing the hydromechanical control of plant growth.

Multicellular organisms grow and acquire their shapes through the differential expansion and deformation of their cells. Recent research has addressed the role of cell and tissue mechanical properties in these processes. In plants, it is believed that growth rate is a function of the mechanical stress exerted on the cell wall, the thin polymeric layer surrounding cells, involving an effective viscosity. Nevertheless, recent studies have questioned this view, suggesting that cell wall elasticity sets the growth rate or that uptake of water is limiting for plant growth. To assess these issues, we developed a microfluidic device to quantify the growth rates, elastic properties and hydraulic conductivity of individual Marchantia polymorpha plants in a controlled environment with a high throughput. We characterized the effect of osmotic treatment and abscisic acid on growth and hydromechanical properties. Overall, the instantaneous growth rate of individuals is correlated with both bulk elastic modulus and hydraulic conductivity. Our results are consistent with a framework in which the growth rate is determined primarily by the elasticity of the wall and its remodelling, and secondarily by hydraulic conductivity. Accordingly, the coupling between the chemistry of the cell wall and the hydromechanics of the cell appears as key to set growth patterns during morphogenesis.

DELLA proteins recruit the Mediator complex subunit MED15 to coactivate transcription in land plants.

DELLA proteins are negative regulators of the gibberellin response pathway in angiosperms, acting as central hubs that interact with hundreds of transcription factors (TFs) and regulators to modulate their activities. While the mechanism of TF sequestration by DELLAs to prevent DNA binding to downstream targets has been extensively documented, the mechanism that allows them to act as coactivators remains to be understood. Here, we demonstrate that DELLAs directly recruit the Mediator complex to specific loci in Arabidopsis, facilitating transcription. This recruitment involves DELLA amino-terminal domain and the conserved MED15 KIX domain. Accordingly, partial loss of MED15 function mainly disrupted processes known to rely on DELLA coactivation capacity, including cytokinin-dependent regulation of meristem function and skotomorphogenic response, gibberellin metabolism feedback, and flavonol production. We have also found that the single DELLA protein in the liverwort Marchantia polymorpha is capable of recruiting MpMED15 subunits, contributing to transcriptional coactivation. The conservation of Mediator-dependent transcriptional coactivation by DELLA between Arabidopsis and Marchantia implies that this mechanism is intrinsic to the emergence of DELLA in the last common ancestor of land plants.

Myosin XI, a model of its conserved role in plant cell tip growth.

In eukaryotic cells, organelle and vesicle transport, positioning, and interactions play crucial roles in cytoplasmic organization and function. These processes are governed by intracellular trafficking mechanisms. At the core of that trafficking, the cytoskeleton and directional transport by motor proteins stand out as its key regulators. Plant cell tip growth is a well-studied example of cytoplasm organization by polarization. This polarization, essential for the cell's function, is driven by the cytoskeleton and its associated motors. This review will focus on myosin XI, a molecular motor critical for vesicle trafficking and polarized plant cell growth. We will center our discussion on recent data from the moss Physcomitrium patens and the liverwort Marchantia polymorpha. The biochemical properties and structure of myosin XI in various plant species are discussed, highlighting functional conservation across species. We further explore this conservation of myosin XI function in the process of vesicle transport in tip-growing cells. Existing evidence indicates that myosin XI actively organizes actin filaments in tip-growing cells by a mechanism based on vesicle clustering at their tips. A hypothetical model is presented to explain the essential function of myosin XI in polarized plant cell growth based on vesicle clustering at the tip. The review also provides insight into the in vivo localization and dynamics of myosin XI, emphasizing its role in cytosolic calcium regulation, which influences the polymerization of F-actin. Lastly, we touch upon the need for additional research to elucidate the regulation of myosin function.

The cAMP signaling module regulates sperm motility in the liverwort Marchantia polymorpha.

Adenosine 3',5'-cyclic monophosphate (cAMP) is a universal signaling molecule that acts as a second messenger in various organisms. It is well established that cAMP plays essential roles across the tree of life, although the function of cAMP in land plants has long been debated. We previously identified the enzyme with both adenylyl cyclase (AC) and cAMP phosphodiesterase (PDE) activity as the cAMP-synthesis/hydrolysis enzyme COMBINED AC with PDE (CAPE) in the liverwort Marchantia polymorpha. CAPE is conserved in streptophytes that reproduce with motile sperm; however, the precise function of CAPE is not yet known. In this study, we demonstrate that the loss of function of CAPE in M. polymorpha led to male infertility due to impaired sperm flagellar motility. We also found that two genes encoding the regulatory subunits of cAMP-dependent protein kinase (PKA-R) were also involved in sperm motility. Based on these findings, it is evident that CAPE and PKA-Rs act as a cAMP signaling module that regulates sperm motility in M. polymorpha. Therefore, our results have shed light on the function of cAMP signaling and sperm motility regulators in land plants. This study suggests that cAMP signaling plays a common role in plant and animal sperm motility.

KATANIN-mediated microtubule severing is required for MTOC organisation and function in Marchantia polymorpha.

Microtubule organising centres (MTOCs) are sites of localised microtubule nucleation in eukaryotic cells. Regulation of microtubule dynamics often involves KATANIN (KTN): a microtubule severing enzyme that cuts microtubules to generate new negative ends, leading to catastrophic depolymerisation. In Arabidopsis thaliana, KTN is required for the organisation of microtubules in the cell cortex, preprophase band, mitotic spindle and phragmoplast. However, as angiosperms lack MTOCs, the role of KTN in MTOC formation has yet to be studied in plants. Two unique MTOCs - the polar organisers - form on opposing sides of the preprophase nucleus in liverworts. Here, we show that KTN-mediated microtubule depolymerisation regulates the number and organisation of polar organisers formed in Marchantia polymorpha. Mpktn mutants that lacked KTN function had supernumerary disorganised polar organisers compared with wild type. This was in addition to defects in the microtubule organisation in the cell cortex, preprophase band, mitotic spindle and phragmoplast. These data are consistent with the hypothesis that KTN-mediated microtubule dynamics are required for the de novo formation of MTOCs, a previously unreported function in plants.

Elucidation of arsenic detoxification mechanism in Marchantia polymorpha: The role of ACR3.

The arsenic-specific ACR3 transporter plays pivotal roles in As detoxification in yeast and a group of ancient tracheophytes, the ferns. Despite putative ACR3 genes being present in the genomes of bryophytes, whether they have the same relevance also in this lineage is currently unknown. In this study, we characterized the MpACR3 gene from the bryophyte Marchantia polymorpha L. through a multiplicity of functional approaches ranging from phylogenetic reconstruction, expression analysis, loss- and gain-of-function as well as genetic complementation with an MpACR3 gene tagged with a fluorescent protein. Genetic complementation demonstrates that MpACR3 plays a pivotal role in As tolerance in M. polymorpha, with loss-of-function Mpacr3 mutants being hypersensitive and MpACR3 overexpressors more tolerant to As. Additionally, MpACR3 activity regulates intracellular As concentration, affects its speciation and controls the levels of intracellular oxidative stress. The MpACR3::3xCitrine appears to localize at the plasma membrane and possibly in other endomembrane systems. Taken together, these results demonstrate the pivotal function of ACR3 detoxification in both sister lineages of land plants, indicating that it was present in the common ancestor to all embryophytes. We propose that Mpacr3 mutants could be used in developing countries as low-cost and low-technology visual bioindicators to detect As pollution in water.

Divergent evolution of the alcohol-forming pathway of wax biosynthesis among bryophytes.

The plant cuticle is a hydrophobic barrier, which seals the epidermal surface of most aboveground organs. While the cuticle biosynthesis of angiosperms has been intensively studied, knowledge about its existence and composition in nonvascular plants is scarce. Here, we identified and characterized homologs of Arabidopsis thaliana fatty acyl-CoA reductase (FAR) ECERIFERUM 4 (AtCER4) and bifunctional wax ester synthase/acyl-CoA:diacylglycerol acyltransferase 1 (AtWSD1) in the liverwort Marchantia polymorpha (MpFAR2 and MpWSD1) and the moss Physcomitrium patens (PpFAR2A, PpFAR2B, and PpWSD1). Although bryophyte harbor similar compound classes as described for angiosperm cuticles, their biosynthesis may not be fully conserved between the bryophytes M. polymorpha and P. patens or between these bryophytes and angiosperms. While PpFAR2A and PpFAR2B contribute to the production of primary alcohols in P. patens, loss of MpFAR2 function does not affect the wax profile of M. polymorpha. By contrast, MpWSD1 acts as the major wax ester-producing enzyme in M. polymorpha, whereas mutations of PpWSD1 do not affect the wax ester levels of P. patens. Our results suggest that the biosynthetic enzymes involved in primary alcohol and wax ester formation in land plants have either evolved multiple times independently or undergone pronounced radiation followed by the formation of lineage-specific toolkits.

Growth control of Marchantia polymorpha gemmae using nonthermal plasma irradiation.

Several studies have documented that treatment by cold atmospheric pressure plasma (CAPP) on plants foster seed germination and growth in recent years. However, the molecular processes that underlie the action of CAPP on the seeds and plants remain mostly enigmatic. We here introduce gemmae of Marchantia polymorpha, a basal liverwort, as a novel model plant material suitable for CAPP research. Treating the gemmae with CAPP for a constant time interval at low power resulted in consistent growth enhancement, while growth inhibition at higher power in a dose-dependent manner. These results distinctly demonstrate that CAPP irradiation can positively and negatively regulate plant growth depending on the plasma intensity of irradiation, offering a suitable experimental system for understanding the molecular mechanisms underlying the action of CAPP in plants.

Unveiling the molecular mechanisms of arsenic tolerance and resilience in the primitive bryophyte Marchantia polymorpha L.

This study addresses the pressing issue of high arsenic (As) contaminations, which poses a severe threat to various life forms in our ecosystem. Despite this prevailing concern, all organisms have developed some techniques to mitigate the toxic effects of As. Certain plants, such as bryophytes, the earliest land plants, exhibit remarkable tolerance to wide range of harsh environmental conditions, due to their inherent competence. In this study, bryophytes collected from West Bengal, India, across varying contamination levels were investigated for their As tolerance capabilities. Assessment of As accumulation potential and antioxidant defense efficiency, including SOD, CAT, APX, GPX etc. revealed Marchantia polymorpha as the most tolerant species. It exhibited highest As accumulation, antioxidative proficiency, and minimal damage. Transcriptomic analysis of M. polymorpha exposed to 40 µM As(III) for 24 and 48 h identified several early responsive differentially expressing genes (DEGs) associated with As tolerance. These includes GSTs, GRXs, Hsp20s, SULTR1;2, ABCC2 etc., indicating a mechanism involving vacuolar sequestration. Interestingly, one As(III) efflux-transporter ACR3, an extrusion pump, known to combat As toxicity was found to be differentially expressed compared to control. The SEM-EDX analysis, further elucidated the operation of As extrusion mechanism, which contributes added As resilience in M. polymorpha. Yeast complementation assay using Δacr3 yeast cells, showed increased tolerance towards As(III), compared to the mutant cells, indicating As tolerant phenotype. Overall, these findings significantly enhance our understanding of As tolerance mechanisms in bryophytes. This can pave the way for the development of genetically engineered plants with heightened As tolerance and the creation of improved plant varieties.

Plant evolution: A tapetum is now effectively present in all land plant lineages.

The tapetum, a tissue that elsewhere ensures correct spore development, is missing in some bryophytes. A new study shows that, in the liverwort, Marchantia polymorpha, a gene controlling spore wall deposition is expressed in the capsule lining, so these cells essentially function as a tapetum.

MiRNAs differentially expressed in vegetative and reproductive organs of Marchantia polymorpha - insights into their expression pattern, gene structures and function.

MicroRNAs regulate gene expression affecting a variety of plant developmental processes. The evolutionary position of Marchantia polymorpha makes it a significant model to understand miRNA-mediated gene regulatory pathways in plants. Previous studies focused on conserved miRNA-target mRNA modules showed their critical role in Marchantia development. Here, we demonstrate that the differential expression of conserved miRNAs among land plants and their targets in selected organs of Marchantia additionally underlines their role in regulating fundamental developmental processes. The main aim of this study was to characterize selected liverwort-specific miRNAs, as there is a limited knowledge on their biogenesis, accumulation, targets, and function in Marchantia. We demonstrate their differential accumulation in vegetative and generative organs. We reveal that all liverwort-specific miRNAs examined are encoded by independent transcriptional units. MpmiR11737a, MpmiR11887 and MpmiR11796, annotated as being encoded within protein-encoding genes, have their own independent transcription start sites. The analysis of selected liverwort-specific miRNAs and their pri-miRNAs often reveal correlation in their levels, suggesting transcriptional regulation. However, MpmiR11796 shows a reverse correlation to its pri-miRNA level, suggesting post-transcriptional regulation. Moreover, we identify novel targets for selected liverwort-specific miRNAs and demonstrate an inverse correlation between their expression and miRNA accumulation. In the case of one miRNA precursor, we provide evidence that it encodes two functional miRNAs with two independent targets. Overall, our research sheds light on liverwort-specific miRNA gene structure, provides new data on their biogenesis and expression regulation. Furthermore, identifying their targets, we hypothesize the potential role of these miRNAs in early land plant development and functioning.

Rapid Propagation of Ca2+ Waves and Electrical Signals in the Liverwort Marchantia polymorpha.

In response to both biotic and abiotic stresses, vascular plants transmit long-distance Ca2+ and electrical signals from localized stress sites to distant tissues through their vasculature. Various models have been proposed for the mechanisms underlying the long-distance signaling, primarily centered around the presence of vascular bundles. We here demonstrate that the non-vascular liverwort Marchantia polymorpha possesses a mechanism for propagating Ca2+ waves and electrical signals in response to wounding. The propagation velocity of these signals was approximately 1-2 mm s-1, equivalent to that observed in vascular plants. Both Ca2+ waves and electrical signals were inhibited by La3+ as well as tetraethylammonium chloride, suggesting the crucial importance of both Ca2+ channel(s) and K+ channel(s) in wound-induced membrane depolarization as well as the subsequent long-distance signal propagation. Simultaneous recordings of Ca2+ and electrical signals indicated a tight coupling between the dynamics of these two signaling modalities. Furthermore, molecular genetic studies revealed that a GLUTAMATE RECEPTOR-LIKE (GLR) channel plays a central role in the propagation of both Ca2+ waves and electrical signals. Conversely, none of the three two-pore channels were implicated in either signal propagation. These findings shed light on the evolutionary conservation of rapid long-distance Ca2+ wave and electrical signal propagation involving GLRs in land plants, even in the absence of vascular tissue.

Convergent evolution of water-conducting cells in Marchantia recruited the ZHOUPI gene promoting cell wall reinforcement and programmed cell death.

A key adaptation of plants to life on land is the formation of water-conducting cells (WCCs) for efficient long-distance water transport. Based on morphological analyses it is thought that WCCs have evolved independently on multiple occasions. For example, WCCs have been lost in all but a few lineages of bryophytes but, strikingly, within the liverworts a derived group, the complex thalloids, has evolved a novel externalized water-conducting tissue composed of reinforced, hollow cells termed pegged rhizoids. Here, we show that pegged rhizoid differentiation in Marchantia polymorpha is controlled by orthologs of the ZHOUPI and ICE bHLH transcription factors required for endosperm cell death in Arabidopsis seeds. By contrast, pegged rhizoid development was not affected by disruption of MpNAC5, the Marchantia ortholog of the VND genes that control WCC formation in flowering plants. We characterize the rapid, genetically controlled programmed cell death process that pegged rhizoids undergo to terminate cellular differentiation and identify a corresponding upregulation of conserved putative plant cell death effector genes. Lastly, we show that ectopic expression of MpZOU1 increases production of pegged rhizoids and enhances drought tolerance. Our results support that pegged rhizoids evolved independently of other WCCs. We suggest that elements of the genetic control of developmental cell death are conserved throughout land plants and that the ZHOUPI/ICE regulatory module has been independently recruited to promote cell wall modification and programmed cell death in liverwort rhizoids and in the endosperm of flowering plant seed.

Transcription of the Antisense Long Non-Coding RNA, SUPPRESSOR OF FEMINIZATION, Represses Expression of the Female-Promoting Gene FEMALE GAMETOPHYTE MYB in the Liverwort Marchantia polymorpha.

Sexual differentiation is a fundamental process in the life cycles of land plants, ensuring successful sexual reproduction and thereby contributing to species diversity and survival. In the dioicous liverwort Marchantia polymorpha, this process is governed by an autosomal sex-differentiation locus comprising FEMALE GAMETOPHYTE MYB (FGMYB), a female-promoting gene, and SUPPRESSOR OF FEMINIZATION (SUF), an antisense strand-encoded long non-coding RNA (lncRNA). SUF is specifically transcribed in male plants and suppresses the expression of FGMYB, leading to male differentiation. However, the molecular mechanisms underlying this process remain elusive. Here, we show that SUF acts through its transcription to suppress FGMYB expression. Transgene complementation analysis using CRISPR/Cas9D10A-based large-deletion mutants identified a genomic region sufficient for the sex differentiation switch function in the FGMYB-SUF locus. Inserting a transcriptional terminator sequence into the SUF-transcribed region resulted in the loss of SUF function and allowed expression of FGMYB in genetically male plants, leading to conversion of the sex phenotype from male to female. Partial deletions of SUF had no obvious impact on its function. Replacement of the FGMYB sequence with that of an unrelated gene did not affect the ability of SUF transcription to suppress sense-strand expression. Taken together, our findings suggest that the process of SUF transcription, rather than the resulting transcripts, is required for controlling sex differentiation in M. polymorpha.

Thermospermine Is an Evolutionarily Ancestral Phytohormone Required for Organ Development and Stress Responses in Marchantia Polymorpha.

Thermospermine suppresses auxin-inducible xylem differentiation, whereas its structural isomer, spermine, is involved in stress responses in angiosperms. The thermospermine synthase, ACAULIS5 (ACL5), is conserved from algae to land plants, but its physiological functions remain elusive in non-vascular plants. Here, we focused on MpACL5, a gene in the liverwort Marchantia polymorpha, that rescued the dwarf phenotype of the acl5 mutant in Arabidopsis. In the Mpacl5 mutants generated by genome editing, severe growth retardation was observed in the vegetative organ, thallus, and the sexual reproductive organ, gametangiophore. The mutant gametangiophores exhibited remarkable morphological defects such as short stalks, fasciation and indeterminate growth. Two gametangiophores fused together, and new gametangiophores were often initiated from the old ones. Furthermore, Mpacl5 showed altered responses to heat and salt stresses. Given the absence of spermine in bryophytes, these results suggest that thermospermine has a dual primordial function in organ development and stress responses in M. polymorpha. The stress response function may have eventually been assigned to spermine during land plant evolution.

Efficient Production of Flavonoid Glucuronides in Escherichia coli Using Flavonoid O-Glucuronosyltransferases Characterized from Marchantia polymorpha.

As a model liverwort, Marchantia polymorpha contains various flavone glucuronides with cardiovascular-promoting effects and anti-inflammatory properties. However, the related glucuronosyltransferases have not yet been reported. In this study, two bifunctional UDP-glucuronic acid/UDP-glucose:flavonoid glucuronosyltransferases/glucosyltransferases, MpUGT742A1 and MpUGT736B1, were identified from M. polymorpha. Extensive enzymatic assays found that MpUGT742A1 and MpUGT736B1 exhibited efficient glucuronidation activity for flavones, flavonols, and flavanones and showed promiscuous regioselectivity at positions 3, 6, 7, 3', and 4'. These enzymes catalyzed the production of a variety of flavonoid glucuronides with medicinal value, including apigenin-7-O-glucuronide and scutellarein-7-O-glucuronide. With the use of MpUGT736B1, apigenin-4'-O-glucuronide and apigenin-7,4'-di-O-glucuronide were prepared by scaled-up enzymatic catalysis and structurally identified by NMR spectroscopy. MpUGT742A1 also displayed glucosyltransferase activity on the 7-OH position of the flavanones using UDP-glucose as the sugar donor. Furthermore, we constructed four recombinant strains by combining the pathway for increasing the UDP-glucuronic acid supply with the two novel UGTs MpUGT742A1 and MpUGT736B1. When apigenin was used as a substrate, the extracellular apigenin-4'-O-glucuronide and apigenin-7,4'-di-O-glucuronide production obtained from the Escherichia coli strain BB2 reached 598 and 81 mg/L, respectively. Our study provides new candidate genes and strategies for the biosynthesis of flavonoid glucuronides.

A subclass II bHLH transcription factor in Marchantia polymorpha gives insight into the ancestral land plant trait of spore formation.

Sporopollenin is often said to be one of the toughest biopolymers known to man. The shift in dormancy cell wall deposition from around the diploid zygotes of charophycean algae to sporopollenin around the haploid spores of land plants essentially imparted onto land plants the gift of passive motility, a key acquisition that contributed to their vast and successful colonization across terrestrial habitats.1,2 A putative transcription factor controlling the land plant mode of sporopollenin deposition is the subclass II bHLHs, which are conserved and novel to land plants, with mutants of genes in angiosperms and mosses divulging roles relating to tapetum degeneration and spore development.3,4,5,6,7 We demonstrate that a subclass II bHLH gene, MpbHLH37, regulates sporopollenin biosynthesis and deposition in the model liverwort Marchantia polymorpha. Mpbhlh37 sporophytes show a striking loss of secondary wall deposits of the capsule wall, the elaters, and the spore exine, all while maintaining spore viability, identifying MpbHLH37 as a master regulator of secondary wall deposits of the sporophyte. Localization of MpbHLH37 to the capsule wall and elaters of the sporophyte directly designates these tissue types as a bona fide tapetum in liverworts, giving support to the notion that the presence of a tapetum is an ancestral land plant trait. Finally, as early land plant spore walls exhibit evidence of tapetal deposition,8,9,10,11,12 a tapetal capsule wall could have provided these plants with a developmental mechanism for sporopollenin deposition.

Conserved CKI1-mediated signaling is required for female germline specification in Marchantia polymorpha.

In land plants, gametes derive from a small number of dedicated haploid cells.1 In angiosperms, one central cell and one egg cell are differentiated in the embryo sac as female gametes for double fertilization, while in non-flowering plants, only one egg cell is generated in the female sexual organ, called the archegonium.2,3 The central cell specification of Arabidopsis thaliana is controlled by the histidine kinase CYTOKININ-INDEPENDENT 1 (CKI1), which is a two-component signaling (TCS) activator sharing downstream regulatory components with the cytokinin signaling pathway.4,5,6,7 Our phylogenetic analysis suggested that CKI1 orthologs broadly exist in land plants. However, the role of CKI1 in non-flowering plants remains unclear. Here, we found that the sole CKI1 ortholog in the liverwort Marchantia polymorpha, MpCKI1, which functions through conserved downstream TCS components, regulates the female germline specification for egg cell development in the archegonium. In M. polymorpha, the archegonium develops three-dimensionally from a single cell accumulating MpBONOBO (MpBNB), a master regulator for germline initiation and differentiation.8 We visualized female germline specification by capturing the distribution pattern of MpBNB in discrete stages of early archegonium development, and found that MpBNB accumulation is restricted to female germline cells. MpCKI1 is required for the proper MpBNB accumulation in the female germline, and is critical for the asymmetric cell divisions that specify the female germline cells. These results suggest that CKI1-mediated TCS originated during early land plant evolution and participates in female germ cell specification in deeply diverged plant lineages.